
THEORY OF NOISE GENERATION FROM 
MOVING BODIES WITH AN APPLICATION 
TO HELICOPTER ROTORS 

F. Farassat 

Langley Research Center 
Hampton, Va. 23663 





NATIONAL AERONAUTICS AND SPACE ADMINISTRATION • WASHINGTON, 0. C. • DECEMBER 1975 


1. Report No. 

NASA TR R-451 


2. Government Accession No. 


3. Recipient's Catalog No. 


4. Title and Subtitle 

THEORY OF NOISE GENERATION FROM MOVING BODIES 
WITH AN APPLICATION TO HELICOPTER ROTORS 


5. Report Date 

December 1975 


6. Performing Organization Code 


7. Author(s) 

F. Farassat, The George Washington University, Joint Institute 
for Acoustics and Flight Sciences 


8. Performing Organization Rep>ort No. 

L-10379 


9. Performing Organization Name and Address 
NASA Langley Research Center 
Hampton, Va. 23665 


10. Work Unit No. 

505-10-26-02 


11. Contract or Grant No. 


12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D.C. 20546 


13. Type of Report and Period Covered 

Technical Report 


14. Sponsoring Agency Code 


15. Supplementary Notes 


16. Abstract 

Several expressions for the determination of the acoustic field of moving bodies are pre- 
sented in this report. The analysis is based on the Ffowcs Williams -Hawkings equation. Apply- 
ing some criteria proposed herein, one of these expressions is singled out for numerical com- 
putation of acoustic pressure signature. The compactness of sources is not assumed and the 
main results are not restricted by the observer position. The distinction between compact and 
noncompact sources on moving surfaces is discussed. Some thickness noise calculations of 
helicopter rotors and comparison with experiments are included which suggest this mechanism 
as the source of high-speed blade slap of rotors. 


17. Key Words (Suggested by Author(s)) 

Aeroacoustics 
Aerodynamic noise 
Helicopter rotor noise 


18. Distribution Statement 

Unclassified - Unlimited 


Subject Category 71 


19. Security Classif. (of this report) 

20. Security Classif. (of this page) 

21. No. of Pages 

22. Price* 

Unclassified 

Unclassified 

61 

$4.25 


For sale by the National Technical Information Service, Springfield, Virginia 22161 


f 


L 


CONTENTS 


Page 

SUMMARY 1 

INTRODUCTION 1 

SYMBOLS 3 

THEORETICAL FORMULATION 8 

General Discussion 8 

Formulation 9 

Solution of Governing Equation 9 

Compact and Noncompact Sources 13 

APPLICATION TO HIGH-SPEED PROPELLERS AND HELICOPTER ROTORS .... 15 

Specialization of General Equation 15 

Thickness Noise of Helicopter Rotors 20 

CONCLUDING REMARKS 22 

APPENDIX - MATHEMATICAL DERIVATIONS AND NOTES 23 

Embedding Procedure 23 

Derivation of Ffowcs Williams -Hawkings Equation by Embedding Procedure ... 27 

Solution of Wave Equation With Generalized Functions as Source Terms 30 

Discussion of Singularities of the Solution of Wave Equation 38 

Relation of S -Surface Solution to Other Forms of the Solution of Wave 

Equation 43 

REFERENCES 46 

FIGURES 48 


THEORY OF NOISE GENERATION FROM MOVING BODIES 


WITH AN APPLICATION TO HELICOPTER ROTORS 

F. Farassat 

The George Washington University 
Joint Institute for Acoustics and Flight Sciences 

SUMMARY 

In this report, several formulations for the determination of the acoustic field of 
bodies of arbitrary shape and motion are presented. The analysis is based on the Ffowcs 
Williams -Hawkings equation which is rederived by an embedding procedure and is given 
in the appendix. The required parameters are the body geometry, time history of motion 
and its surface pressure distribution. Based on criteria proposed in the report, one of 
these formulations which consists of an observer time derivative and integrals over the 
curve of intersection of a collapsing sphere and the body, is singled out for numerical cal- 
culation of acoustic pressure signature. Compact or noncompact source distributions on 
a moving body are defined in terms of the time scale of fluctuations of the sources and the 
length scale of a virtual radiating surface whose geometry depends on the body shape and 
motion. A new expression is derived for the computation of the acoustic pressure signa- 
ture of rotating blades. This expression is not restricted by compactness or far -field 
assumptions. A careful study of singularities of the various forms of the solution is 
undertaken and the nature of these singularities is explained. Some examples of thickness 
noise of helicopter rotors in hover and forward flight are worked out. These examples 
suggest strongly that high-speed blade slap of rotors is mainly the thickness noise 
phenomenon. 


INTRODUCTION 

The generation of sound from moving bodies, such as a rotating propeller or heli- 
copter rotor blades, is an undesirable byproduct of many useful machines. The estimation 
of the sound field produced by the motion of such bodies is, therefore, essential in design- 
ing quiet machines. The theory for estimation of this sound has been developing gradually 
in recent years. Generally, one is concerned with the motion of a slender body such as a 
fan or a propeller blade with a rather complicated flow pattern around it. The pressure 
disturbance of the fluid medium caused by the motion of the body is felt as sound when the 


disturbance passes over any observer in the medium. However, without proper approxi- 
mations, one is faced with a difficult theoretical problem in calculating the acoustic 
pressure. 

An important approximation is introduced by assuming that the body does not disturb 
the medium appreciably and thus the nonlinearities may have negligible effects. The gov- 
erning equation for the acoustic pressure is then a wave equation with some inhomogene- 
ous source terms depending on the net local force of the body on the fluid and the normal 
velocity of the surface of the body. Historically, the effects of these two terms were 
studied separately. 

Gutin (ref. 1) was first to calculate the sound from static propellers. He employed 
a result derived by Lamb (ref. 2) for the acoustic field of a stationary concentrated force. 
Because of the periodicity of the propeller forces, one is able to replace them by stationary 
and harmonically oscillating forces which are distributed on the propeller plane. Gutin 
applied Lamb’s result to these forces. For harmonic analysis the governing wave equa- 
tion was converted into the Helmholtz equation whose solution gave the acoustic pressure 
spectrum. 

Garrick and Watkins (ref. 3) removed Gutin' s restriction of static propellers by 
first studying the acoustic field of a concentrated force in uniform rectilinear motion. 

The observer is in the frame fixed to the source. In the nonrotating frame moving with 
the propeller, the analysis of Gutin was repeated and the acoustic pressure spectrum was 
obtained. 

Deming (ref. 4) studied the effect of the source term involving the normal velocity 
distribution of the body surface for static propellers. He utilized Rayleigh's relation for 
the piston radiation into a semi-infinite space and obtained the acoustic pressure spectrum. 
This noise is referred to as the thickness noise. Arnoldi (ref. 5) extended Deming's work 
to propellers in uniform axial motion. Lyon (ref. 6) has studied the thickness noise of 
helicopter rotors in forward flight. 

One of the common assumptions in treating the problem of acoustic radiation is the 
compactness of the sources. A stationary source is compact if its dimension is much 
smaller than the wave length of radiation. This definition needs to be modified for moving 
sources. Effectively a compact source may be treated as a point source and considerable 
simplification is introduced in the acoustic analysis (ref. 7). Even the forward flight 
effects of helicopter rotors may be studied numerically rather easily in this case (ref. 8). 

It is the main purpose of this report to develop a theory for the calculation of the 
acoustic pressure signature for bodies in arbitrary motion in the ground-fixed coordinate 
system where the observer is not limited to the far field and no compactness assumption 
is made. The starting point of the present analysis is the equation developed by Ffowcs 
Williams and Hawkings, denoted the FW-H equation, using the acoustic analogy (ref. 9). 
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With a view toward application to helicopter rotors and propellers, the second-order 
source term which depends on Lighthill's stress tensor is neglected. The resulting equa- 
tion is valid in the entire unbounded space and contains source terms which involve Dirac 
delta functions. These terms indeed come from the boundary conditions of the problem. 
One may then solve the FW-H equation by using the Green's function for the unbounded 
space. 

An interesting fact about the FW-H equation is that it is possible to obtain many dif- 
ferent forms of the solution (ref. 9). All these forms have some type of singularities for 
bodies moving at high speeds. In this report several forms of the solution of the FW-H 
equation are presented. The criteria which are proposed later will be applied to single 
out one form for the computation of the acoustic pressure signature. The distinction be- 
tween the compact and noncompact source models for moving bodies is discussed in some 
detail. In many situations of interest, one cannot assume that the sources that generate 
the sound are compact. 

A new expression for the determination of the acoustic field of high-speed helicopter 
rotors and propellers is derived and presented here. The part of this expression called 
the thickness noise is used for some trend calculations in the case of helicopter rotors. 
These calculations indicate that the thickness noise may be the mechanism of high-speed 
blade slap of helicopters. Some comparison with experiments is given to support this 
proposition. 

The mathematical derivation and notes presented in the appendix are supplementary 
to the main body of the report. The technique of embedding a problem with a restricted 
domain into a problem with unbounded domain is discussed first. This technique is em- 
ployed to derive the FW-H equation. The different forms of the solution of the FW-H 
equation are obtained and their relation to each other is explained in the appendix. A 
special emphasis is put on the problem of singularities in the solution of the FW-H equa- 
tion, The nature of these singularities is also discussed in the appendix. 


SYMBOLS 


c 


speed of sound 


D 


see equation (30) 


DP 


disk plane of helicopter rotor on propeller 



viscous stress tensor 
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F(y; x,t) = f(y,t-r/c) = 0 equation of S -surface 


f( ) arbitrary function 

f(y = 0 equation of surface of body 

G(x,^),G(x,|) Green's functions 


g 

= T - t + r/c 

H() 

Heaviside function 

h( ) 

arbitrary function 


camber surface of blade in -frame 

L(x,t) 

length scale of 2 -surface for a fixed observer position x and time t 

1 

length scale of the body 

Mn 

Mach number based on velocity normal to body, v^/c 

Mr 

Mach number based on velocity in radiation direction, Vr/c 

Mt 

tip Mach number 

m -curves 

system of curves on 2 -surface orthogonal to F -curves 

Ni 

unit outward normal to 2 -surface 

"i 

unit outward normal to body 

PiJ 

compressive stress tensor 

P,Pl,P2 

pressure, acoustic pressure 

pfp 

pressure distribution on body due to thickness 
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Q(y ,r),Qj(y ,t) Source distribution functions 


^i 


net force of body on fluid 


r.n 


distance between source point and observer 


ret 


subscript denoting retarded time 


unit vector in radiation direction 


components of unit vector f in ?f' -frame 


^min 


minimum distance from observer to 2 -surface 


see figure A2 


S 

S* 


surface of body f(y ,t) = 0 

part of body surface whose sources contribute to acoustic pressure at position 
X and observer time t 


Tij 


Lighthill stress tensor 


T(y ) arbitrary bounded function 


thickness distribution function of blade in rf' -frame 


observer time 


fp>W 


time scales of pressure fluctuations and variation of speed of body 


u; 


fluid velocity 


un 


fluid velocity normal to the body 


vehicle velocity 
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V[ component of vehicle velocity V in -frame 

n = -V’l + 

= -v^ - 

Vj velocity of body 

Vi(?,T) velocity of a point of the body f(y ,t) =0 in Lagrangian coordinates 

Vn component of velocity normal to body surface, v • n 

Vj. component of velocity of body surface in radiation direction, v • r 

Xj observer position vector 

y -frame Cartesian coordinate system fixed to the undisturbed medium 
yj source position vector 

r curve of intersection of body and sphere g = 0 

Ap pressure distribution on rotor blade due to lift 

6( ) Dirac delta function 

6ij Kronecker delta 

? Lagrangian coordinate of a point on the body 

t],?]' Cartesian coordinate frames 

6 angle between n^ and r^ 

6 angle between Ni and the radiation direction rj 

angle between normal to camber surface and rj 
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I 2 

= (l + Mn'^ - 2Mn cos Q) 
density 

density of undisturbed medium 

= P -Po 

surface generated by F(y; x,t) = 0 
source times 

T*,T* emission time of source 

T(j(x,t) time spent for sphere g = 0 to cross body for fixed observer position x 

and observer time t 

tq see figure A2 

(p angle between tjj- and ?? 2 -axes 

0(x,t),0 j(x ,t),02(^>t) unknown functions in wave equations (see appendix) 


A 

P 

Po 

P 

s 


SI 


angular velocity of rotor 


□ 


1 o ^ o 

wave operator, -2— - 

C^ 




Laplacian operator, 


8^ + 9^ + _9 
9yi^ 9 y 2^ 9ys 


Subscripts 1, 2, and 3 are used to denote components of a vector or partial deriva- 
tives along the corresponding axis. 

Vectors are also denoted by an arrow over the symbol. 

Summation convention on subscripts i,j is used throughout this report, 
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The element of the volume in each reference frame is written in terms of the differ- 
ential of the symbol used to designate the frame (i.e., dy is the element of volume in the 
y -frame). 


THEORETICAL FORMULATION 


General Discussion 

For a body in motion, whose surface is described as f(y,T) = 0, the governing equa- 
tion for the determination of the acoustic pressure p is the Ffowcs Williams -Hawkings 
(FW-H) equation (ref. 9): 


^ 6 ( 0 ] 



6(f) 

9yj 


9yi9yj 


( 1 ) 


The summation convention on subscripts i and j is used here. The derivation of this 
equation is based on the acoustic analogy and it is, therefore, valid in the entire unbounded 
space. The y -frame in which equation (1) holds is fixed with respect to the undisturbed 
medium. The above equation is derived in the appendix. 

Equation (1) is the generalization of the equation derived by Lighthill (ref. 10) in 
connection with his jet noise theory. The first two inhomogeneous source terms in equa- 
tion (1) are absent in Lighthill's theory and the contribution from the remaining term 
a^Tjjyayjayj is the so-called quadrupole noise due to turbulence. The first term in the 
above equation arises as a result of the motion of the surface in the normal direction. 

Each surface element can be viewed as a small piston acting on the fluid with speed v^. 
The second term comes from the local surface stress which consists of viscous 

stress Ejj and the thermodynamic pressure p. Physically, this term is the contribu- 
tion of the net force acting on the fluid due to the viscous stress and the pressure distribu- 
tion on the body surface. 

It is assumed that the flow parameters on the moving surface are known. These 
parameters are obtained experimentally or by aerodynamic calculations. Obviously, this 
assumption simplifies the acoustic analysis considerably. As in the case of jet noise the- 
ory, not all the features of the acoustic field are expected to be exhibited by the present 
formulation. This is because of the fact that in the situations of interest, the flow around 
the body is quite complex and the required flow parameters could be specified to a limited 
degree. Despite this shortcoming, the acoustic analogy appears as one of the most useful 
methods for the study of the noise from bodies in motion. 
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Formulation 


Since in the cases of interest, one is concerned with the motion of a slender body 


such as a rotor or fan blade, the region of turbulent flow is small and of relatively low in- 
tensity. Therefore, the term involving Tij in equation (1) will be neglected. It is well 
known that turbulence is a very inefficient noise producing mechanism (ref. 10). This is 
yet another reason to neglect the last term of equation (1) which is a volume source. Tur- 
bulence has another effect on the acoustic field by producing fluctuating pressure on the 
body surface which is taken care of analytically through the second term in equation (1). 
The viscosity effect will also be neglected. Under the above assumptions, the governing 
equation for the generation of the sound for moving bodies is: 


Here, n is the unit outward normal to the body. The pressure on the right side of the 
equality is the surface pressure. In the solution of this equation, the geometry and the 
time history of the motion of the body will be assumed known. By the nature of the acous- 
tic analogy, the surface pressure distribution is a given parameter of the problem. 


In view of the fact that it is possible to get numerous forms of the solution of the 
FW-H equation, some criteria for the selection of a suitable form for the computational 
work are needed. The following criteria are therefore proposed: 

(a) The solution should require as little information as possible about the source 
distribution on the surface of the body. 

(b) The singularities in the solution should be integrable. 

(c) The computation time should be as low as feasible. 

The solution of equation (2) is p = Pj + P 2 where pj and P 2 are the solutions 
of the following equations: 



|^[poVnlvf| 6(f)] -^[pnilVfl 6(f)] 


( 2 ) 


Solution of Governing Equation 



(3) 



(4) 
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These equations are of the types discussed in the appendix (eqs, (A31) and (A32), respec- 
tively). In the following paragraphs the solutions of equations (3) and (4) are written with 
the above criteria in mind. 

Consider the surface F(y ; x ,t) = f(y ,t-r/c) = 0. Here, x and t are the observer 
position and time, respectively. These two parameters are kept fixed in the following dis- 
cussion. It is shown in the appendix that the surface F = 0, denoted by S, serves as a 
virtual radiating surface for the body f(y,T) = 0 which is in motion. The 2 -surface is 
the locus of the points on the body whose signals arrive simultaneously to the observer at 
time t. This surface is generated by the curves of the intersection of the body and the 
collapsing sphere g = t - t + r/c = 0 for -oo < t ^ t. These curves of intersection are 
called the F-curves here. The collapsing sphere in the process of crossing a two-bladed 
rotor system is shown in figure 1. (See also fig. (Al).) In terms of the S -surface, the 
solution of equation (3), using equation (A34) with Q(y,r) = PqVq, is 


47t pi(x,t) = 



d2 


(5) 


/ 2 

The symbol A is used for (1 + - 2Mfj cos dj where = v^/c and 9 is the 

angle between n and t (=x - y , y is the source position). The subscript ret 
stands for retarded time. The singularities in the integrand of equation (5), which appear 
when A = 0, are shown to be integrable in the appendix. The reason for keeping the time 
derivative outside the integral rather than using the result similar to equation (A44) is 
twofold. First, the order of singularity in the integrand is increased when the derivative 
is taken inside the integral. Second, as seen from equation (A44), one is required to cal- 
culate more parameters on the S -surface than those needed in equation (5). Also, when 
equation (A44) is expanded by performing the differentiations in the integrand, the result 
is a very long and complicated expression which does not seem to be useful for numerical 
computation. 

For equation (4), the solution equivalent to equation (A46) with Q^Cy ,t) = pnj, is 


4tt P2(x,t) 



( 6 ) 


In the discussion following equation (A46), the reason for preferring this form of solution 
to those involving space derivatives, such as equation (A36), is stated. The remarks con- 
cerning equation (5) apply equally here. 
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The solution of equation (2), from equations (5) and (6), is 


4tt p(x,t) = 47 t(pi + pg) = 1 ^ 


r l /PpCVn + P cos 
Jf= 0 ^ \ A /ret 



(V 


Although the 2 -surface is defined anal 3 d;ically by F(y; x,t) = f(y ,t-r/c) = 0, in gen- 
eral the above equation cannot be integrated in closed form even if the surface pressure 
p is known analytically. This is because of the implicit representation of the Z -surface 
which may have a complicated geometrical shape in the situations of interest. (See, for 
instance, figs. A3 and A5.) Because of this difficulty one needs to introduce a change of 
variables of integration to get other forms of the solution. However, one must be aware 
of the additional singularities introduced in the solution by such transformations. When 
the relation of equation (A68) is employed, the following equivalent forms of equation (7) 
are obtained: 


, a C^2 r PnCVn + P cos 0 „ C^2 C cp cot 9 

Att p(x ,t) = \ ^ \ ^ dr dr + — dr dT 

at Jjj Jp r sin 0 Jp ^.2 


C at Jo* \ r 1 - Mr ^ "^5 


p cos 9 


dS 


S* 1 - M, 


( 8 ) 

(9) 


As described earlier, r in equation (8) is the curve of intersection of the collaps- 
ing sphere g = t - t + r/c = 0 and the body f(y, t) = 0. This sphere enters and leaves 
the body at the source times and T 2 . If, for a fixed x and t, the collapsing 
sphere enters and leaves the body more than once, then the integration in equation (8) is 
carried out over all the periods for which the intersection occurs. 

Equation (9) applies to rigid bodies only. The time r* is the emission time of the 
sources on the element of surface area dS of the body. The region of the body surface 
which contributes to the acoustic pressure at x and at the time t is denoted by S*. If 
all the points on the body surface move at speeds less than the speed of sound c, then S* 
is the entire surface of the body. Otherwise, one has S* = S*(5T,t). In the preceding equa- 
tion, Mj. is defined as v^/c where v^ = v • f and v is the local velocity of the body 
surface. The time derivatives in equations (7), (8), and (9) should be taken numerically in 
the applications. 

For the application of equation (8), one needs to construct the F-curves. This is, 
however, relatively simple. At any source time Tj, one starts with the largest sphere 
with its center at the observer position x which touches the body. Let the radius of this 
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sphere shrink at the speed of sound. The body should continue its motion. The intersec- 
tion of the sphere and the body forms the T-curves, The sphere is allowed to collapse 
until it leaves the body completely at the time Tg. The observer time corresponding to 
this sweep is t = rj + rj/c where rj is the radius of the sphere at the source time 
T 2 . One must make sure that the sphere does not intersect the body prior to Tj and 
after t 2 . This may be figured by knowing the time history of the motion of the body. 
Otherwise, the contribution of further intersections to the integrals in equation (8) must 
be included in the calculation of the acoustic pressure. The singularities of equation (8) 
are discussed in the appendix. They are integrable and could be dealt with in numerical 
work. 

If a Lagrangian coordinate system ^ is introduced for the points on the body, the 
following transformation rule exists between the y - and ^ -frames: 

y = f f v(f ,r) dT = y (f,T) (10) 

^0 

Here it is assumed that the two frames coincide at t = 0. The velocity of a point on the 
body with coordinate ^ and at the time t is v(^ ,t). The emission time t* is the 
solution or the solutions of the equation 


g= T* 



( 11 ) 


for fixed x , t, and ? . In practice, this equation is transcendental and rarely can a 
simple analytic solution be obtained. However, the above equation can usually be solved 
numerically by an iteration technique. It must be mentioned that, formally, one may write 

T* = t - = t - 

c c 


Because of this relation, it is customary to use the retarded time notation for formulas 
with integrals over the body surface. Such a notation could lead to confusion when taking 
a derivative inside an integral. For this reason, to convert the space derivatives in equa- 
tion (A36) to a time derivative, the L -surface solution is used in the appendix even if the 
solution consisting of an integral over the body surface could also be employed. In the 
latter case the conversion is not as straight forward as that based on the S -surface 
solution. 

For numerical calculations based on equation (9), the body surface S is first 
divided into a number of small elements. For each element, equation (11) should be em- 


12 


ployed to evaluate the corresponding source time t* (keeping x and t fixed) at an 
appropriate point on this element. However, because t* needs to be calculated a large 
number of times for each t, the use of equation (9) may require an excessive amount of 
computation time. The problem of the singularities of this equation is not easy to handle 
in numerical work. 

Going back to the criteria proposed earlier in this section, it is seen that the most 
suitable form of the solution of the FW-H equation, equation (2), for computational pur- 
poses is given by equation (8). This solution satisfies criterion (a) by requiring only the 
surface pressure and not the pressure gradient. The choice of F and r as the vari- 
ables of integration is to satisfy the criteria (b) and (c). Also only the observer time de- 
rivative and not the derivatives with respect to Xj (such as in eqs. (A3 5) and (A36)) 
appears in equation (8). This reduces the computation time and thus favorably affects 
criterion (c). 


Compact and Noncompact Sources 

The terms compact and noncompact sources are used very frequently in acoustics. 
The purpose of this section is to bring out more clearly the distinction between compact 
and noncompact sources on moving bodies. Using the S -surface solution, one can define 
these terms precisely. 

To be specific, the following equation will be considered: 

tfp = -^[pntlvfl 6(f)] (12) 

This equation is similar to equation (A32) with Qj(y ,t) = -pnjlvfl 6(f). The solution given 
by equation (A36) is 


47t p(x ,t) = dS 

Jf=0 ' ^ ^ret 


(13) 


After using the relation of equation (A68), one obtains 


4n p(x,t) = 



(14) 


where r* is the solution of equation (11) for points on the surface S of the body. 

For a fixed observer position x and time t, let the S -surface have the largest 
dimension L(x ,t) and let the time that the sphere g = 0 takes to cross the entire body 
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be T(>(x ,t). If the time scale tp of the pressure variation on the surface S and the 
time scale ty of the variation of speed of the body are much greater than T(.(x ,t) and 
if the shortest distance from the observer to the S -surface r^in is much greater than 
L(x,t), then equation (14) can be approximated as follows: 


47t p(x ,t) = 



dS = - 



(15) 


where q is the net force on the medium due to the pressure on the surface S and t* 
is the mean source time which is obtained by solving 




- t 


+ 


X 




= 0 


(16) 


Here the mean location of the body is given by yg(r*). It is obvious that, effectively, the 
body is treated as a point source and in fact equation (15) is the solution of 



which is what Lowson solved and published in 1965 (ref. 7). 

It is seen that the crucial assumptions in treating the body as a point source are 

^p ”^c(^ >1) 

Tv » Tc(x ,t) V (18) 

^^min L(x ,t)^ 

Under these conditions, the source distribution is called compact. Otherwise, it is called 
noncompact. These conditions are dependent on observer time and position. Therefore, 
if an observer time interval is considered, all the above conditions should hold for this 
interval. Because of the dependence on observer position, the term "compact source dis- 
tribution" is misleading since compactness is not a characteristic of the source distribu- 
tion per se. For low subsonic speeds of the body, the length scale L(x ,t) and the typical 
length of the body I are of the same order of magnitude so that the third condition of 
equation (18) can be replaced by rjjjin 
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The most important fact about compact source representation is that one is able to 
relate the acoustic field to global parameters of the motion of the body such as the net 
force q of the body on the fluid or in the case of thickness noise to net rate of mass in- 
jection. However, in the case where the source distribution is noncompact, the local pres- 
sure on the S -surface at the source time t (=t - r/c) and not the net force on the body 
at the time t contributes to the acoustic pressure at the observer. 

For rotating bodies, even at moderate tip speeds, some of the conditions of equa- 
tions (18) could not be satisfied for the observer in or close to the plane of rotation. This 
may explain the poor agreement of the compact source theories with experiments so far. 

APPLICATION TO HIGH-SPEED PROPELLERS AND HELICOPTER ROTORS 

Specialization of General Equation 

The solution of equation (2) is used here to derive a useful expression for high-speed 
propellers and helicopter rotors. The particular solution employed is equation (8). The 
main assumption here is that the intersection of the collapsing sphere g = 0 with the 
blade can be approximated by the curve of intersection of g = 0 with the disk plane. Be- 
cause of this approximation, it is suggested that the observer be located in the region be- 
tween 45° above or below the plane of rotation during the period when g = 0 crosses the 
system. This is the region where the noncompact source formulation should be used so 
that this requirement does not affect the usefulness of the results. The observer position 
is otherwise arbitrary. It can be in either the near or the far field. 

A single blade is considered first. Three coordinate systems are used as shown in 
figure 2. The first which is called the -frame is fixed to the blade such that the V\n2~ 
plane coincides with the plane of rotation. The 772 -axis is along the span and the rj'^-axis 
is parallel to the chord of the blade. The origin of this frame is at the rotation center. 

The second system is the t] -frame with its origin at the same location as the -frame 
but is nonrotating. The 772772 -plane coincides with the disk plane. The axes of this frame 
are parallel to those of the third system, the y -frame, which is fixed to the undisturbed 
medium. The position vector of the origin of the rf -frame and the If' -frame is 4 (r) and 
the angle between the 77*2 -axis and the 77 -axis is <p. In the following discussion, the com- 
ponents of any vector in the if' -frame are primed for distinction. The vehicle speed with 
respect to the y -frame which is d| /dr is denoted by V. 

The coordinate transformations between the preceding frames are 

if = y - F (19) 
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77*1 = cos 0+772 


772 = -77j sin 0 + 772 cos 0 > 

^3 = ^3 

v\ = (yi - 4i) cos 0 + (y 2 - ^ 2 ) si” 

^2 ^ "(yi ■ ^ 1 ) 4> + {Y2 - ^ 2 ) 

^3 = y3 - ^3 


( 20 ) 


( 21 ) 


Let 773 = ^{VifV 2 ) and ^?3 = ’^ip\fV2) equations of the camber surface and 

the thickness distribution, respectively. The equation of the surface of the blade in the 
-frame is f{rf') = 0 where 


i(v') 


% - ^(^ 1 .^ 2 ) - 

V^3 - '^('?i>'?2) + 



( 22 ) 


Here the abbreviations US and LS stand for the upper surface and the lower surface of the 
blade, respectively. The preceding equations are written in such a way that f (if') < 0 
inside the blade and f (Jf') > 0 outside the blade. The equation of the blade surface in the 
y -frame is 


f(y,T) = f [jT(y,T)] = 0 


(23) 


where if'{y,T) is given by equations (21). 

The normal velocity v^ of the blade surface is found from the relation 

-8f/3T 


( 24 ) 
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From equation (23) one can obtain the following results: 


9t 


'(Ti + hi)(V’i - 77^0) + (T2 + h2)(V^ + V\^) - 
Fl - hi)(V’i - v'2^) + (T2 - h2)(V^ + 77’ifi) + 


(US)'l 

S> (25) 

(LS) 


|vf|2 = 


l + (Ti + hi)2 + (T2 + h2)^ 


1 + (Ti - hi)2 + (T2 - h2)^ 




(US) 

(LS) 


(26) 


where are the components of the vehicle velocity V in the ?f' -frame, SI is 

the angular velocity of the blade system, and the subscripts 1 and 2 on the functions T 
and h designate partial derivatives with respect to rj'-^ and respectively. Since 
the blade may be considered a slender body, one can assume that 


|hi| « 1 

iTil « 1 


|h2| « 1 |T2| « 1 


(Except near leading edge) 


(27) 


For this reason, one obtains 

|vf|^ = 1 + 

over the upper and lower surfaces of the blade. 

Equation (23) also may be employed to obtain the relation 


Vf • r = |vf| cos 9 = 


f-(Ti + hi)f'i - (T2 + h2)r^ + 


L-(Ti - hl)?i - (T2 - h2)?^ - 


(US)' 

(LS) 


(28) 


(29) 


where (r'i,r 2 >r 3 ) s^re the components of the unit radiation vector f in the if' -frame. 
Combining the preceding equation and equations (27) and (28) yields 


|vf|^ sin2 0 = 1 - r^2 ^^2^^ _ j.^2) = £,2 


(30) 
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Equation (8) will be rewritten as follows: 


p(x ,t) = I: dr dT + 4 f P - 9°^ ^ dr dr+ r 9E-£^drdT 

9t Jtj Jr r sm 0 9t Jxj Jr Jtj Jr 


(31) 


For the term involving Vq, which is in fact the solution of equation (3), equations (24), 
(25), and (30) give 


Pl(x,t) 


^J.r2r Vn dr - f '^1^1 ^ '^2^2 

4 n atJjj Jr r sin 0 2 ir 9 t Jj^^ -^r(DP) J'D 


dr dT 


(32) 


where Vj and V 2 are given by the following relations: 
V’l = -V\ + t?2^ 

> 

= -V^ - 


(33) 


The abbreviation DP stands for the disk plane of the helicopter rotor or propeller. 

For the terms which involve p on the right side of equation (31), consider the 
following integral: 


I = r P_92il dF dT 

Jri Jr r 

When equations (29) and (30) are used, one obtains 

I = r (Pu-PL)(-hin -^2^2-^ ^ 3 ) dr dr 

Jtj Jr(DP) rD 


^■^2 r ( Pu + Pl)('^1^1 '^2^^2 ) 
Jti Jr(DP) rD 


(35) 


where the subscripts U and L stand for the upper and lower surfaces of the blade, re- 
spectively. Using the notation Ap h pj^ " Pu Pt ^ (Pu ''' Pl)/^’ equation (35) is 

written as 
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(36) 


I 


pr dr dT - 2 P f dr dr 

'^Tj ^r(DP) rD Jri Jr(DP) 


where 0^ is the angle between the normal to the camber surface (pointing into the region 
7/3 > 0) and the radiation direction. In equation (36), the parameters Ap and p,p are 
pressure distribution due to lift and thickness, respectively. Since, in the situation of in- 
terest, these two parameters are of the same order of magnitude, the two integrals in 
equation (36) are of the same order of magnitude if the observer is in or very near the 
77 i? 72 -plane when the sphere g = 0 crosses the blade. Directly above and below this 
plane, the first integral in equation (36) dominates. 

The contribution from the terms containing pressure on the right side of equa- 
tion (31) is, therefore. 


2 at ^r(DP) 4?7rD 


pr 

-’r(DP) 


Ap cos 0h + ^ 2 ^ 2 ) 

47rr^D 


dP dT 


(37) 


This is in fact the solution of equation (4). For several blades the acoustic pressure from 
each blade must be added linearly in equations (32) and (37). Essentially a similar result 
consisting of partial derivatives with respect to observer position rather than time t was 
obtained and reported in reference 11 . 

The contribution to the acoustic pressure denoted by pj(x ,t) will be called the 
thickness noise although for a blade with zero lift, Ap = 0, the acoustic pressure has a 
nonzero contribution from P 2 since p-j, 0 on the blade. This definition is consistent 
with what has been called the thickness noise in the literature. (See refs. 4, 5, and 12.) 

There is some evidence (ref. 6 ) that in or close to the 77 ^ 772 -plane, the plane of rota- 
tion, the contribution of thickness noise pj is greater than P 2 at high tip speeds. This 
fact and the fact that the calculation of pj only requires knowledge of blade geometry and 
kinematics led to the development of a computer program to calculate the thickness noise 
(ref. 14). At present, the lack of sufficient information about the surface pressure distri- 
bution on rotating blade systems prevents a reliable calculation of P 2 . However, the out- 
put of the thickness noise program and a comparison with experimental results may be 
employed to find under what conditions the acoustic pressure p can be approximated by 

Pi- 
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Thickness Noise of Helicopter Rotors 

To study the thickness noise of high-speed helicopter rotors, equation (32) was used 
to develop a computer program for a rotor system with rectangular blades and a uniform 
airfoil section along the span. The program and the user's guide are presented in 
reference 14. 

In this program, the F-curves are approximated by the arcs of intersection of the 
collapsing sphere g = 0 and the disk plane which is in the r]^ri 2 -pla.ne. The double inte- 
gral of equation (32) is approximated by a finite difference technique with a smaller mesh 
size near the leading edge of the blade. The derivative with respect to the observer time 
is taken numerically after the result of the integration is smoothed to reduce enlargement 
of errors. As expected, due to the construction of many F-curves for each observer time, 
the computation time is very long. A special effort should be made to design an efficient 
algorithm for the computation of noise from noncompact source distributions. 

The computer program was used to investigate the relation between thickness noise 
and high-speed blade slap of helicopter rotors. The particular features of blade slap of 
concern here are 

(a) Impulsive acoustic pressure signature 

(b) Strong dependence of sound intensity on tip Mach number 

(c) Directionality in the plane of rotation 

The following numerical examples are presented to point out that thickness noise 
indeed exhibits all the above characteristics of high-speed blade slap. The rotor system 
for figures 3 to 5 has two blades with biconvex parabolic arc airfoil sections. The rotor 
diameter is 10 m and the blade chord is 0.4 m. The thickness ratio of each blade is 
10 percent. The observer is in the rotation plane 50 m from the center of rotation. For 
figures 6 and 7 the rotor parameters are those of a test helicopter. The two-bladed rotor 
diameter is 13.42 m and the blade chord is 0.68 m. The thickness ratio of the blades is 
9.3 percent. The rotor speed is 324 rpm. The helicopter is at an altitude of 15 m and the 
observer is 80 m ahead of the helicopter. 

Figure 3 shows the theoretical variation of the acoustic pressure signature of a 
hovering helicopter as the tip speed is increased. The observer is in the plane of rotation. 
As the tip Mach number is increased, one sees a steepening of the signature, becoming im- 
pulsive above Mj; = 0.6. The shape of the acoustic pressure signature changes at and 
above sonic tip speed with the new shape having two positive peaks and a negative peak of 
approximately equal magnitude. Between M^ = 0.6 and near sonic speed, the negative 
peak is dominant. The absolute value of the negative peak varies as a function of tip Mach 
number according to exp(8.94Mt^) for this particular configuration in the range 


20 


Mt = 0.4 to 1.1. Therefore, both features (a) and (b) of the blade slap are exhibited by the 
thickness noise. 

Figure 4 shows the theoretical variation of the acoustic pressure signature with the 
observer elevation. The observer distance is 50 m from the center of rotation of the 
hovering rotor system. The tip Mach number is kept fixed at Mt = 0.9. The direction- 
ality of the signal in the plane of rotation and the reduction in the impulsiveness of the 
signature at higher elevation is observed in this figure. Feature (c) of high-speed blade 
slap is, therefore, displayed by the thickness noise. 

Figure 5 shows that it is the region near the tip of the rotor that generates most of 
the noise. The cut-out of the blade system of figures 3 and 4 is increased from 14 per- 
cent to 84 percent. Comparing this figure with figure 3(f), one can see that the pressure 
signature has not changed significantly either in shape or level. 

Figure 6 shows the theoretical effect of change in the thickness distribution on the 
acoustic pressure signature. Note that in this figure, the rotor system is different from 
previous figures and a steady forward speed is also included. It is seen that the changes 
in the pressure signatures are significant and should be considered in designing a low- 
noise profile for the region near the blade tips. From equation (32), it may be shown that 
the acoustic pressure for the same section type (e.g., NACA four -digit airfoil or super- 
critical airfoil) is proportional to the thickness ratio. Although the supercritical airfoil 
section under the assumed conditions generates a higher noise level, one may use blades 
of smaller thickness ratio, as compared with other types, and obtain good aerodynamic 
characteristics as well as low noise level. 

Figure 7 shows comparison of theoretical and experimental pressure signatures for 
the rotor system of figure 6 with an NACA four -digit airfoil section. The theoretical 
pressure signature is due to thickness alone but the agreement, particularly at higher for- 
ward speed, is very good. The most probable source of discrepancy at lower speeds is 
the neglect of the pressure sources on the blade surfaces of the rotor. 

Since thickness noise does indeed display all the observed features of high-speed 
blade slap of helicopters, including the correct level of acoustic pressure, the thickness 
noise is proposed as the mechanism of the slap phenomenon. The theory and the numeri- 
cal examples suggest several methods of controlling this noise: 

(1) Reduce the tip speed of the rotor as much as possible (see fig. 3) 

(2) Modify the tip region by 

(a) Reducing the thickness ratio; the acoustic pressure is directly proportional 

to thickness ratio (see eq. (32)) 

(b) Changing the thickness distribution of the airfoil section (see fig. 6) 

(c) Altering the planform 
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The alteration of the planform is not included in the computer program presented in ref- 
erence 14. However, its effect on the thickness noise can be suspected by the changes in 
the 2 -surface it produces for different observer times. This effect may be studied 
numerically. 


CONCLUDING REMARKS 

Several expressions for the determination of the acoustic pressure signature of 
bodies in motion are presented herein. Although all these expressions are of theoretical 
value, in practical situations, one form, consisting of an observer time derivative and in- 
tegrals over the curve of intersection of a collapsing sphere and the body, was found to be 
the most useful for numerical work. No assumptions regarding the compactness of the 
sources or the location of the observer were made. An expression for the calculation of 
the acoustic pressure signature of rotating blades was derived. Using this expression, 
some thickness noise calculations were presented which, after comparison with experi- 
mental results, indicated that the thickness noise should be suspected as the cause of the 
high-speed blade slap of helicopter rotors. Further experiments are needed to verify 
this result. 

Lar^ley Research Center 

National Aeronautics and Space Administration 
Hampton, Va. 23665 
October 31, 1975 
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MATHEMATICAL DERIVATIONS AND NCXTES 
Embedding Procedure 

As mentioned in the Introduction, to use the Green's function for the unbounded 
space, the original problem, restricted to the domain outside of the moving body, is em- 
bedded into a class of problems defined in the entire unbounded space. The merits of this 
technique will be discussed shortly. The procedure will be illustrated by an example 
using an ordinary differential equation. 

Consider the differential equation 

y" - y = f(x) 
with boundary conditions 

y(a) = Ya 

y 0 as X ^ oo 

This problem can be solved by finding the solutions of the following equations: 

y'i - Yi = f(x) 

with boundary conditions 
yi(a) = 0 

yj - 0 as X — 00 


/x in [a,oo) 

\f(x) — 0 as X — oo 

> (A 2 ) 


/x in [a,oo) \ 

\f(x) — 0 as X — 00/ 


(Al) 


y*2 - Y2 = 0 

with boundary conditions 
Y2(a) = Ya 
y2 ^ 0 as x — oo 





> (A 3 ) 


J 


The solution to the original problem is then y = y^ + Y2 due to linearity of equations (Al). 
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The Green's function for equations (A2) is (ref. 15) 

G(x,4) = sinh (a - x) H(^-x) + sinh (a - 0 H(x-^) 

where H(x) is the Heaviside function defined by 


H(x) = 


(x > 1)| 

(x < 0)J 


The solution for y is then the following: 
y(x) = (G(x,^),f(^)> + y 2 (x) 

G(x,^) f(|) d^ + y 2 (x) 

= e^”^ C sinh (a - f(^) d^ + sinh (a - x) C e^'^ f(^) d^ + y^ e^"^ 

'^a. *^x 



(A4) 


(A5) 


The solution of equations (A3) which is 


V2 = ya 


(A6) 


has been utilized in equation (A5), 

For use in the upcoming discussion the following relation is obtained by evaluating 
the derivative of the solution of y(x) in equation (A5) at x = a: 


ya'^yk^'f e^"^ f(|) d| 

^a 

where y^ = y’(a). Now consider the following differential equation: 

y" - y = h(x) 

with boundary conditions 

y — 0 as |x| — 00 


fx in (- 00 , 00 ) \ 

lh(x) - 0 as |x| - 00 / 


(A7) 


(A8) 
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The Green's function for this problem is 


G(x,^) = H(^x) + e-''+^ H(x-^)] 


The solution of equation (A8) is, therefore, 


y(x) = (G(x,^),h(|)> 



(A9) 


Let y(x) be the solution of equations (Al) corresponding to a given function f(x) 
satisfying the condition f(x) — 0 as x — oo. The function y(x) = H(x-a) y(x) certainly 
satisfies the boundary conditions of the equations (A8). What is the corresponding source 
function h(x)? This function turns out to be a generalized function. It is obtained as 
follows: 


y"(x) - y(x) = H(x-a) [y"(x) - y(x)] + 6(x-a) y^ + 6'(x-a) ya 
= H(x-a) f(x) + 6(x-a) y^ + 6'(x-a) ya 

where 6(x-a) is the Dirac delta function. And, thus, the function h(x) is given by 


y'(x) = H(x-a) y’(x) + H’(x-a) y(x) 


= H(x-a) y'(x) + 6(x-a) y(x) 


= H(x-a) y'(x) + 6(x-a) Ya 


y"(x) = H(x-a) y"(x) + 6(x-a) y'(x) + 6'(x-a) ya 


H(x-a) y"(x) + 6(x-a) y^ + 6'(x-a) ya 


h(x) = H(x-a) f(x) + 6(x-a) y^ + 6'(x-a) Ya 


(AlO) 
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The function h(x) is not determined completely since Ya is known but ya is not given 
as seen from equations (Al). Assume that y^ is either measured by an experiment or 
is otherwise known. It will be calculated shortly. Note that the domain of the definition 
of h(x) is the entire unbounded real line. What has been achieved is that the problem 
described by equations (Al) may be viewed as a special case of the class of problems de- 
scribed by equations (A8). This technique which extends the domain of the definition of a 
differential equation will be called the embedding procedure. It is clear that if the Green's 
function for the problem with the larger domain is known, then the new problem may be 
solved assuming that the appropriate additional boundary conditions, such as y^ in equa- 
tion (AlO), are known. Since in many problems, particularly those involving partial dif- 
ferential equations, it is hard to find the Green's function, the embedding procedure ex- 
tends the usefulness of the available Green's functions. 

To convince the reader that one obtains the same solution of equation (A5) for the 
problem of equations (Al), the source function h(x) in equation (AlO) is substituted in 
equation (A5). Some elementary properties of the generalized function theory (ref. 16, 
pp. 1-25) will be used here. One obtains 


Since one is interested in the value of y when x is in [a,oo), from equations (Al) the 
following is obtained: 


H(x-a) y(x) = y(x) 


= <G(x,|),h(|)> 

= <G(x,^),H(^a) m + 6(1 -a) y^ + 6'(|-a) y^) 


= (G(x,^),H(|-a) f(^)> + y^ G(x,a) - y^ 



e^"^ H(x-a) + H(x-a) 


(All) 



a g-x+a 


(A12) 
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Letting x = a, precisely the same relation as equation (A7) is obtained. This indicates 
that given y^, one cannot choose y^ arbitrarily, but it is obtained uniquely in the pro- 
cess of solving the problem. The relation between the boundary conditions, equation (A7), 
will be called the compatibility condition. Substituting y^ in terms of y^ from equa- 
tion (A7) in equation (A12), one obtains the same solution as equation (A5). 

Looking back, it is seen that the solution of the problem described by equations (Al) 
with restricted domain was obtained by first embedding the problem in the class of prob- 
lems described by equations (A8) with unbounded domain. This provided the inhomogene- 
ous source term which is the generalized function h(x) involving conditions at the point 
X = a not entirely supplied by the original problem, equations (Al). However, the method 
eventually led to the compatibility condition (A7) which yielded the unique solution to 
equations (Al). 

For the Ffowcs Williams -Hawkings (FW-H) equation, the compatibility condition is 
in fact an integral equation in terms of the surface pressure of the body which is usually 
difficult to solve. One may then use experimentally obtained surface pressure distribu- 
tions. The usefulness of the method presented here should not be overlooked as it pro- 
vides a systematic and neat method of attack which may work. 

Derivation of Ffowcs Williams -Hawkings Equation by Embedding Procedure 

Let the surface of the moving body be given by the function f(y,r) = 0. The func- 
tion f(y ,t) is defined to be positive outside and negative inside the body. Following the 
procedure of the preceding section, all the flow parameters will be multiplied by the func- 
tion H[f(y, t)], where H(f) is the Heaviside function, in order to embed the problem of 
the external flow and acoustics of this body into the class of problems valid in the entire 
unbounded space. 

It is required on physical grounds that all flow parameters vanish at infinity. For 
this reason, the quantity P - Pq where Pq is the density of the undisturbed medium, 
will be used as the density parameter, rather than p itself. However, if p is multi- 
plied by a quantity that vanishes at infinity, such as the flow velocity, it is kept unchanged. 
Outside the body, the continuity and momentum equations hold true. Using summation 
convention on subscripts i and j, these equations are 

Continuity equation: 

^(P - Po) + = 0 (A13) 

or dy| 
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Momentum equation: 


■* ^ '"I")) 


= 0 


(A14) 


where uj is the fluid velocity, t is the time, and P^j is the compressive stress ten- 
sor. The parameter p - Pq will be designated p. The frame of reference y is fixed 
with respect to the undisturbed medium. 

The continuity equation, valid in the entire unbounded space, is 


Br 


[p H(f)] 




H(f)] = H(f) 


Bt 8yj^ 


^\bt ^ 8yJ 8yj 



+ u 


8H 


+ PqU^ — 

° ' 9Yi 


(A15) 


Note that 






= 6(f) ^(ui - Vi) 

= 6(f) |vf|(un - Vn) = 0 


(A16) 


where vi is the velocity of the body and the no-penetration condition un = Vn is used 
on the body. The subscript n denotes the component normal to the body. Using the pre- 
ceding result, equation (A15) will be written in the following form: 

lf[2 ««)] * Wi 


= PoUn|vf| 6(f) 


= Po^nlvfl 6(f) 


(A17) 
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Using the above technique, the momentum equation can be written as: 


_ 8 _ 

ar 


[pui H(t)] + + pulUj) H(t) 



= PlJ ^ e«) (A18) 

In this relation, the result of equation (A16) has been applied. Taking the time deriva- 
tive of both sides of equation (A 17) and the space derivative with respect to of both 
sides of equation (A 18) and then subtracting the latter from the former, one obtains 


^[p H(f)] - 


9yi9yjL 


(Pij + puiuj) H(f)j = |-[pQVn|vf| 6(f)] - ^ 


’ij If. 


(A19) 


Now add and subtract the term c^ a^j^jo H(f)Jy/ayi2 on the left side of equation (A19) and 
rearrange the term involving Py to get 



c2-9?^[p H(f)] =|-[poVn|vf| 6(f)] 
9yi^ 




9yi9yj 


(Pij + puiUj - c^p 6ijj H(f) 


(A20) 


Here c is the speed of sound in the undisturbed medium. The quantity 

Tjj = Pjj + pujuj - c^p 6ij (A21) 

is called the Lighthill stress tensor, well known in jet noise theory. Let p = c^p. This 
quantity is the acoustic pressure in the region not too close to the body. It is the pertur- 
bation pressure in linear theory even close to the body itself. Using this notation, equa- 
tion (A20) will be written as 

□2[p H(f)] = - v2j[p H(f)] 


= |^[poVn|vf| 6(f) 


a 

9yi 




9yi9yj 


Tij H(f) 


(A22) 
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The above equation was derived first by Ffowcs Williams and Hawkings (ref. 9). 

They based their derivation of the generalized mass continuity equation (eq. (A17)) and 
momentum equation (eq. (A18)) on the study of mass and momentum flow through a con- 
trol volume, respectively. The body is represented by a surface of discontinuity in the 
unbounded medium. Inside the body, the fluid is at rest with the conditions of the undis- 
turbed medium. The method presented in this appendix was known to Ffowcs Williams 
and Hawkings and the technique, without using the terminology of embedding, was sketched 
in their original work (ref. 9). Note that the FW-H equation is valid even if shock discon- 
tinuities are present in the flow field. This is because of the fact that this equation is 
based on conservation of mass and momentum which should hold in the presence of shock 
waves. 

It should be emphasized that as in the case of ordinary differential equations, the 
source terms in equation (A22) can not all be specified arbitrarily. For example, the sur- 
face pressure on the body should be specified in such a way that if the observer is located 
on the surface of the body at the time t, the calculated acoustic pressure p should be 
the same as the local surface pressure at the time. This is the compatibility condition. 

When obtaining the solution of equation (A22), the function which is solved for will 
be written as p(x,t) rather than p(x,t) H(f). The latter function is actually what the 
solution represents. Mathematically, the solution is valid anywhere including inside the 
body for which p(x ,t) H(f) = 0. This fact can be used to advantage in some problems. 
When the formal solution of equation (A22) is regarded as an integral equation in terms of 
the body surface pressure, one can show that H(f) = 1/2 on the surface of the body. 

Solution of Wave Equation With Generalized Functions 
as Source Terms 

Let f(y ,t) =0 be the equation of the surface of a moving body. The fu ;tion 
f(y,r) is defined such that its gradient Vf points outside the body. The soh ion of the 
inhomogeneous wave equation 

2 

1 _ v20 = Q(y ,T) |vf| 6(f) (A23) 

will be considered first. Next the solutions of the wave equation with the following source 
terms 


|p[Q(y,r) Ivfl 6(f)],^[Qi(y,r) |vf| 6(f)] 
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will be given. In the preceding equation, the functions Q(y,r) and Qi(y,r) are given 
piecewise smooth functions on the body surface. Note that the choice of the preceding 
source terms corresponds to those of FW-H equation (A22). 

The Green's function for the solution of the wave equation in the unbounded space is 
6(g) /47rr where g = t - t + r/c, t is the source time, t is the observer time, 
r = |x - y I and x and y are the observer and source locations respectively. When 
X and t are fixed and the source time t varies from -oo to t the surface g = 0 
is a sphere with center at x whose radius collapses at the speed of sound. The solution 
of the equation (A23) is thus 

4tt 0(x,t) = |vf[ 6(f) 6(g) dy dr (A24) 

The volume integral extends over the unbounded three-dimensional space and the time in- 
tegral is over the interval -OO T = t. 

The integrals involving a product of two delta functions are studied in the literature 
of generalized functions (ref. 16, pp. 209-247). The result is an integral over the subspace 
where f = 0 and g = 0. To get this integral, the variables of integration in the preced- 
ing equation should be transformed to a coordinate system which has f and g as new 
variables. The details are outlined in reference 17. One writes: 


_ df dg dyi 

9(f,g)/9(y2»y3) 


df dg dy^ 
(Vf X vg) . ei 


c df dg dr 

|vf| sin 0 


(A25) 


where ej is the unit vector along yj^-axis and F is the curve of intersection of the sur- 
faces f = Constant and g = Constant. The angle 0 is between Vf and F (=x - y). 

In the preceding equation one may use other combinations of axes to get the same final re- 
sult. Substituting equation (A25) in equation (A24) and integrating with respect to f and 
g, one obtains 


4tt 


<p(x,t) = { ^ ^ 




cQ(y,r) 


r r sin 0 


dr dT 


(A26) 
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Here, because of the products of delta functions in the integrand of equation (A24), T is 
the curve of intersection of the surfaces f = 0 and g = 0. (See figs. 1 and Al.) The 
source times and T 2 correspond to the times when the sphere g = 0 enters and 
leaves the body, respectively. The condition 0 = 0 will be discussed separately later 
on. 

Another form of the solution of equation (A23) may be written by noting that the sub- 
space for which f = g = 0, keeping x and t fixed, is the surface 

F(y; x,t) = f(y,t-g = 0 


This equation is obtained by solving for t from g = 0 and then substituting in 
^(y = 0* This surface, designated as S, is the locus of all the points on the surface 

f(y = 0 whose signals arrive simultaneously at the observer location x and at the 

time t. To visualize this surface, consider a large sphere enclosing the body whose 
radius collapses at the speed of sound to radius zero at observer time t with the center 
at observer position x. As the sphere contracts, let the body trace its history in time. 
The surface generated by the curves of intersection of the sphere and the body is in fact 
the 2 -surface. Incidentally, these curves of intersection are the P -curves employed in 
the equation (A26). Note that for points on a P-curve, or a point lying on the S -surface, 
with coordinate y , the source time t is related to t by the relation ^ = t - ^ where 
r = |x - y |. In figure Al the situation is illustrated for the moving body f(y ,t) = 0. The 
sphere and the body at the source time t are shown with dotted lines. Their intersec- 
tion, the P-curve corresponding to the time t, lies on the 2 -surface. This surface is 
formed by all the P-curves for -oo < t s t. The 2 -surface is not necessarily simply 
connected. 

Now that the 2 -surface has been defined, use the relationship 6(F) 6(g) = 6(f) 6(g) 
in the integrand of equation (A24). This follows from the fact that the subspace defined 
by f = 0 and g = 0 is the same as one defined by F = 0 and g = 0, once again keep- 
ing X and t fixed. Equation (A24) may be written as 

47T 0(x,t) = |vf| 6(F) 6(g) dy dr (A27) 

The domain of the volume integral is the unbounded three-dimensional space and that of 
the time integral is -oo < t ^ t. The following result is employed to integrate the above: 
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d d^ , dyi dy; dF dg 
3(F,g)/8(y3,T) 

_ dyj dy2 dF dg 
aF/9ys 

_ dy^ dy2 dF dg 
N3IVFI 


Here d2 is the element of the surface of F(y, x,t) = 0 and N is the outward unit nor- 
mal to this surface. The symbol 9(F,g)y^9(y3,r^ stands for Jacobian. 

The following expression is obtained from the relation 


F(y; x,t) = f(y,t-g = = 0 


where the subscript ret stands for retarded time: 


cE dF dg 
|VF| 


(A28) 



ret 


In this ejqjression, Mn = Vn/c, Vn = -(9f/9T)/|vf | is the normal velocity of the body, and 
6 is the angle between the outward normal to the body and the vector in the radiation 
direction r . The symbol A will be used in this report as 



(A29) 


When equation (A28) is substituted in equation (A27) and the result integrated with respect 
to F and g, one obtains a new form of the solution for equation (A23). This form is 
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where the integral is over the S -surface (x and t are kept fixed). This result was 
originally published in reference 18. It was derived from equation (A26) using geometric 
arguments. 

From equation (A30), it is seen that the S -surface serves as a virtual radiating body 
whose surface structure depends on the body shape and motion as well as the observer 
position and time. For very small values of M^, the 2 -surface may be approximated by ■ 

the surface of the body itself. From the definition of A, one can show that the condition 
A = 0 occurs only when Mn = 1 and 6 = 0, simultaneously. This is equivalent to the < 

body having sonic speed at the point of tangency with the sphere g = 0. It will be shown 
later that the singularity A = 0 is integrable in the situations of interest. 

Next the solutions of the following equations will be studied; 



(A31) 



(A3 2) 


Since the operator commutes with operators 8/8r and 9/9yi, from equa- 

tions (A26) and (A30), it follows that 



(A33) 



(A34) 



(A3 5) 



(A36) 


Equations (A33) and (A3 5) may be written in other forms by taking the derivatives 
inside the integrals. The method is presented in reference 17. However, the resulting 
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equations may not be suitable for numerical integration in the case of complex geometries. 
The derivatives in equations (A34) and (A36) may also be taken inside the integrals. The 
technique will only be worked out in detail for equation (A34) as follows: 


- 

at Jf=o ^ 


Q(y,T) 


A 


dZ 


ret 


at Jf=o 


Q(y;x,t) dS 
^ IvFl 




' C ft (r ) * C 7 f (A37) 


where Q(y ; x,t) = |Q(y,T) (vf[j . In the preceding equation the following identity, in 


'ret 

which T(y ) is an arbitrary function, is used: 


r T(y) dZ = |VF| T(y) 6(F) dy 

Jp— Q J -oo 


(A38) 


To get this useful identity, note that 


_ dF dyi dy2 
aF/ays 


dF dZ 
|VF| 


(A3 9) 


Substitute this relation on the right side of equation (A38) and then integrate with respect 
to variable F. 

Applying equation (A38), one integral in equation (A37) may be written as 




6(F) dy = f i 

Jf=0 ^ 


aQ 

at 


dZ 

VF| 


(A40) 


The following identity, closely related to one used by Ffowcs Williams and Hawkings 
(ref, 9, eq. 6.2), is required for the second integral on the right side of equation (A37): 
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T(y) 6 ’(F) dy 


■f — 

8F T(y) 

dE 

'^F=o ^yi 

|vf |2 

|VF| 


(A41) 


The summation convention on i is used in equation (A41). Derivation of this identity is 


not trivial. 


It is obtained by writing 



where N is the unit normal vec- 


i=l 

tor to the S -surface. The symbol dy^ is defined as the element of the area of the coor- 
dinate plane whose normal is parallel to yj-axis, that is, dy^ = dy 2 dy 3 , etc. Substituting 
this relation in equation (A39) and then using the result on the left side of equation (A41), 
one obtains 



T(y) 6'(F) dy 


y T(y ) Nj 

J-oo A, |VF| 
1=1 


6'(F) dF dyi 


^+00 
\ > 


T(y ) Ni 

J-oo ^ 9F 

i=l 

Yi 

_ |VF| _ 


6(F) dF dyi 


(A42) 


where 

9F 


indicates that while taking the derivative 8/aF, the coordinates in the plane 

yi 


in which dyi is defined must be kept fixed. For example, 

or 

dinates y 2 and y 3 are kept fixed. It may be shown that 


means that the coor- 


yi 


a 

8F 


f. 9F/9yi 9yi |VF| Ni 9yi 


(i = 1, 2, 3; no sum on i) 


and also dyi/Ni = dS (i = 1, 2, 3; no sum on i). When these two results are used in 
equation (A42) and the latter integrated with respect to F, the result is equation (A41). 

The last integral of equation (A37), by the use of the identity in equation (A41), 

yields 


r q^6 

J-oo r 8t 


(F)dyi=-C Jl/iLJF 
‘ ^F=0 9yi \9yi 9t 



(A43) 
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Equation (A34) can thus be written as 




c 

1 3Q 

-§-i 

(dF 

Q N-^ 

dZ 

^F=0 

r at 

9yi' 

lat 

r|vF| y 

|VF| 


(A44) 


Here the relation Ni = (9F/3yi)/|VF| is used. 

When precisely the same technique is applied to equation (A36), one obtains 


4ir 


02(x,t) = r 

I 


F=0 



_a_( 

aF 

dZ 

9yi' 

y9Xj r|vF| y 

|VF| 


(A45) 


where Qj = (|vf|Qj)^^^. 

The preceding equation will be written in yet another form which is very useful in 
applications. First, note that 


a (Qi) 

Qr 

1 9Qr 

axiVr / 

■ '72 

' cr at 


and 


3F _ ^i 8F 

axj ' c at 


where Q^. = Q^r^ When these are substituted in equation (A45), the following is obtained: 


4tt 


02(x,t) = -f ^ 

Jf=o r2 


dS 

IVFI 



1^ 

r at 


N-1 

9yi\9t r|vF| 7 


dS 

|VF| 


Now compare the second integral on the right side of this equation with the integral in 
equation (A44). Since Q is arbitrary in equation (A44), one may write 


C 

1 9Qr 

--^1 

faF 


dZ 

Jf=o 

r at 

9yi 

i^at 

r|vF| 7 

|VF| 


± r Qr dZ 
at Jjr_Q r lyFl 
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Note that |vF| = (|^f|^)ret' When the above relation is substituted in the equation pre- 
ceding it, the new form of solution for equation (A36) becomes 


47t = 


r ^ dz 

'^F=0 r2 


VF 


1±C Qr dZ 
c at Jf=o ^ |vf| 



(A46) 


where Qj. = Qi?!, 

The importance of the last result lies in the fact that now only one derivative, that 
is, the time derivative, appears in the solution which may be performed numerically for 
a fixed observer position. This is considerably simpler than using equation (A36) in 
which the space derivatives must be calculated numerically at a fixed observer time. 

Equation (A46) may also be written as 


4. 02(5f,t) = - 


^2 r cQr 


r^ sin 9 


dF dr + 


atJ^- Jp 


Q 


^ dr dr 


r sin 9 


(A47) 


To get this from equation (A46), compare equation (A26) with equation (A30). It can also 
be obtained by using equation (A68). 

It is seen that one may write the solution of the wave equations considered in this 
section of the appendix in many different forms. The solution of the FW-H equation sat- 
isfying the criteria proposed in this report is based on equations (A33) and (A47). 

Discussion of Singularities of the Solution of Wave Equation 

One disturbing feature of the different forms of the solution of the wave equation is 
the appearance of singularities such as A = 0 in equation (A34). These singularities 
limit the applicability of the solution. It may appear that such singularities are genuine. 
In this section, it will be shown that, provided that the right variables are chosen in the 
solution of equation (A23), the singularities are integrable. In the case of the following 
expressions or other equivalent forms: 



(A48a) 
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h 




d2 


(A48b) 


it is not permissible to take the derivatives inside the integral without proper interpreta- 
tion, The order of the singularity in the integrand increases when the derivatives are 
taken inside. One should then use Hadamard's method of selecting the finite part of a 
divergent integral or its modern equivalent form in generalized function theory. Due to 
multidimensionality of the integrals, the use of such a technique is an unnecessary com- 
plication which could be avoided by applying numerical differentiation with respect to the 
observer time. The ensuing discussion will, therefore, be confined to the study of the 
singularities of the following integral: 



(A49) 


In this section, it will be shown that the singularities of equation (A49) are integrable. 

Since in the following discussion x and t are kept fixed, the dependence on these 
variables is suppressed in the notation and the integrand will be written as 
T(y) = Q(y ; x,t)/r. Assume that T(y) is bounded on the S-surface, Thus the following 
equivalent form of equation (A49) will be considered for the study of the singularities: 



since VF = |vf|A j.g^ and Q= |vf|Q The set of points on the 2 -surface where 
|vF| =0 is called the critical points of the surface. At these points, provided that f = 0 
has no critical points, = 1 and 0 = 0. This means that the surfaces of f = 0 and 
g = 0 are tangent, and normal body velocity is equal to the speed of sound c. In the fol- 
lowing discussion, it is assumed that not all the second partial derivatives with respect to 
variables y^ are zero at critical points of F(y; x,t) = 0. The following three cases are 
distinguished. 

Case (i). - If the sphere g = 0 is tangent to the body for an extended period of time, 
when the body is on one side of the sphere, the 2 -surface corresponding to this period de- 
generates into a curve in space. Since this curve has no surface content, the integral of 
equation (A50) is zero. This curve itself can become an isolated point in space if g = 0 
and f = 0 become tangent and immediately separate. Again one obtains 1 = 0 in 
equation (A 50), 
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Case (ii).- There may be isolated points on the S -surface where |vF| = 0. This 
situation can occur, for example, at points where the sphere g = 0 enters or leaves the 
body at the speed of sound. Such points may also appear when g = 0 is in the process 
of crossing the body. To visualize this latter condition, consider a finite wing with a rec- 
tangular planform and a symmetric section moving in a straight path with variable speed. 
(See fig. A2(a).) Let the observer be in the plane of symmetry of the wing as shown in 
figure A2(b). If the body accelerates to sonic speed, and subsequently decelerates with- 
out spending some time at sonic speed, the S -surface corresponding to t = tq + Tq/c 
has a singularity at point B. This surface is shown in figure A3. 

This type of sir^ularity appears most often in calculations for rotating bodies. In 
such a case, the tangent plane to the S -surface at the point B is not defined. In fact, in 
the vicinity of B, the S -surface looks like a cone since it can be represented by the func- 
tion F(y ) = Fij-y^yj =0 in a coordinate system with origin at the point B where 
F[j = 8^F /^y^ayjlg. For a small patch of area near B on 2 -surface, designated as D, 
one has 


f T(y ) ^ = T(C ) f 

^ Jd |vf| Jd Ivf| 

for some ^ in region D. Thus, consider the following integral near point B: 

Let the yjy 2 -plane coincide with the tangent plane to the sphere g = 0 when it is 
at the point B. In this plane, set a polar coordinate (R,^) such that 


yj = R cos 


(A 52a) 


y 2 = R sin <p 


Then equation (A51) becomes 
j, ^ r RdRdcp 

|9F/9y3| 


(A52b) 


(A53) 
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and from F(y) = F^jy^yj = 0, one obtains 


9F 


^ = FisYi + F23y2 + 

9Y3 


(A 54) 


Solving for Y3 from the equation of the surface near B, and substituting in equa- 
tion (A54), and finallY using the relations of equations (A52), the following result is 
obtained: 


pRo r 
Jn Jt 


dcp dR 


^ \J2A cos^ cp + B sin 2<p + 2C sin^ cp 


= Ro L 


dcp 


'(p \j(A + C) + (A - C) cos 2cp + B sin 2cp 


(A55) 


where 


A = (fi3 - FiiF33)/2 
B = (F13F23 - F[2F33)/2 



The domain of <p, denoted bY F)(p, is the set of values of cp which makes the denomina- 
tor in the integrand of equation (A55) real. This domain is not necessarilY 2ir. Even 
when the denominator in equation (A55) is zero, there is not a genuine singularitY since it 
behaves as \J\<p - <Pq\ where (Pq is the value of cp which makes the denominator zero. 
If the region D is taken in such a waY that its projection in the YiY2-Pl3^*^® in a circle 
of radius Rq, one obtains I' Rq, Therefore, bY taking Rq small enough, one can 
neglect a patch of area near the tYpe of singularitY discussed here without introducing 
large errors in calculations. 

Case (iii) .- If in the processes of crossing the bodY, the sphere g = 0 becomes tan- 
gent to the bodY at a point whose normal velocitY is sonic and remains sonic for extended 
periods of time, the S -surface intersects itself along a curve. This situation is uncom- 
mon for rotatir^ bodies. Let the speed varY as shown in figure A4, and using the wing of 
figure A2, let the sphere be tangent to the wing at the point B when t = Tq. The 2- 
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surface is shown in figure A5. The structure of this surface near point B is presented 
in figure A6. 

Again, a new y -frame with origin at point B will be introduced with the yiY 2 - 
plane coinciding with the tangent plane to the sphere g = 0 when it is at point B. This 
plane is tangent to the P-curves which cross each other passing through point B. The 
slopes iri]^ and m 2 of the tangents to P-curves in yiy 2 -plane are given by the solution 
of the equation 


F22m2 + 2Fi2m + Fii = 0 


(A56) 


where 


TT. . — 9 




9yi9yj 


B 


The analysis of case (ii) can be applied on that part of the 2 -surface where y 3 < 0. (See 
fig. A6.) This shows that, for this region of the surface, P exists and can be made arbi- 
trarily small if the region D is taken as small. 

For the intersecting part of the S -surface in region y 3 > 0, one can set up a new 
y -frame at any point P on curve BC with the yiy 2 -plane coinciding with the tangent 
plane to the sphere g=0 at P. The normals to the two intersecting parts of the £ - 
surface are defined and could be obtained by using equation (A56) with Fij evaluated at 
point P. One obtains the following relations on each intersecting part of S -surface, say, 
in A’BC near point P 

|VF| = y . vIvFlIp + 0(yi2) (A57a) 

vIvFlIp = Ni VFi|p (A57b) 

where Fj = dF/dy^ and Nj is the unit normal to A'BC. Using equations (A57), one 
can write 


r 


= lim 
e-0 



dF dy3 
QfP + /3y3 


(A58) 
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where a and /3 are related to the second partial derivatives of F(y ; x,t) at point P 
and 6 > € is a small fixed parameter. The above integral exists and can be made as 
small as one wishes by taking the value of the integration limit, 6 small. This means 
that one may neglect a strip of surface along the intersecting part of the 2 -surface which 
includes the curve of intersection BC without producing large errors in calculation. 


The F-curves and the system of curves on the S -surface orthc^onal to them, denoted 
by m -curves, form a grid on the 2 -surface. Except at points on this surface where 0 = 0, 
one has the relation 


At points where 0 = 0, the F-curves or m -curves intersect themselves as shown in fig- 
ure A7. Figure A7(a) illustrates the case when g = 0 enters or leaves the body. Fig- 
ure A7(b) is for the condition 0=0 when the sphere g = 0 is in the process of crossing 
the body. The third case (fig. A7(c)) corresponds to figure A6, when the 2 -surface inter- 
sects itself. In any case, one may neglect a small patch of area near the points for which 
0 = 0, even when Mn = 1. 

In figure A8 a region of the 2 -surface, where 9*0 and * 1, is shown. It is 
seen that 


Relation of 2 -Surface Solution to Other Forms 


of the Solution of Wave Equation 


d2 = dF dm 


(A59) 


dm = -£^ 


(A60) 


sin 0 


where 0 is the angle between Nj and the radiation direction rj. But 



and 



A 


(A61) 
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thus 



(A62) 


Therefore, from equations (A59) to (A62), one gets 



(A63) 


and 


^ _ cA dP dT 


sin 6 


(A64) 


To get a relation between dZ and the element of area of the body surface dS, con- 
sider a portion of the S -surface corresponding to the observer time t. (See fig. A8.) It 
is assumed that the body is rigid. Point A which lies on the P -curve at source time t, 
denoted as r(r), would move to point B on the body at source time t + dr. Note that 
V is the local velocity of the body at the point A when r(r) was formed. Let C be 
the point of intersection of the m -curve from A and the curve r(r + dr). Let F and 
m be the unit tangents to F-curve and m -curve at C, respectively. Then dS is given 


by 


dS = I CB X df 


(A65) 


where dF = F dF. It can be seen that 


CB = -dm + V dr 


Here dm = m dm. Substituting for CB into equation (A65), one gets 



(A66) 


From figure A9 and the definition of dZ, the following may be shown 


-m X F(dm dF) = N dZ 


= nl^MnLdZ 


(A67) 


A 
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In the last step of equation (A67), N is substituted from equation (A61). Also since 
r = r X n/sin 0, it may be shown that 


V X r(dP dr) = V X (r X 

sin 6 

= (vn? - Vrn 

= (Mnr - Mrn)^ (A68) 

where = v • rf, Vr = v • f , and Mr = Vr/c, Here equation (A64) is used to give 

dP dr _ dZ 
sin 6 cA 


From equations (A65), (A67), and (A68), the relation between dZ and dS is ob- 
tained which is 


dS = |l - Mrl ^ 
I r| ^ 


(A69) 


For use in different forms of the solution of the wave equation, equations (A64) and 
(A69) give 


^ = dS (^70) 

A sin 0 |i _ Mr| 

It is seen that when Mr = 1, one gets a singularity in the solution in the form of an 
integral over the body surface S. This condition occurs when |v | = c. When this hap- 
pens, the r -curves corresponding to the same observer time t intersect one another or 
become tangent. These curves may cover a portion of the surface S more than once. 
For example, if the body is moving at steady rectilinear supersonic speed, then only that 
portion of the body which is within the forward Mach cone from the observer will be cov- 
ered twice by the F -curves. Note that only for the solution involving an integral over the 
body surface, it is assumed that the body be rigid. The other forms of the solution apply 
for general surfaces which deform while in motion. 
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Figure 1.- Formation of a F-curve by the intersection of the collapsing sphere 
g = 0 and the body (a helicopter rotor). 


^2 



Figure 2.- Coordinate frames used in determination of the acoustic field 

of rotors and propellers. 
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Pressure 



(a) Mt = 0.4; period, 113.6 msec. 



(c) Mt = 0.6; period, 76.0 msec. 
Figure 3. 




(d) Mt = 0.7; period, 65.1 msec. 

- Theoretical variation of the acoustic pressure signature of a hovering helicopter. 
Thickness noise only. (Note change of scale.) 
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(g) Mt = 1.0; period, 45.5 msec. 


Figure 3. 




Concluded. 


Pressure, N/m 



Figure 4.- Theoretical variation of the acoustic pressure signature with the observer elevation for a hovering 
helicopter. Period, 50.6 msec. Thickness noise only. (Note the change of scale.) 




Blade 


Figure 5.- Influence of the region near the blade tip on the generation 
of thickness noise of a hovering helicopter. Period, 50.6 msec. 
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(b) NACA four -digit airfoil. 



(a) Biconvex airfoil (parabolic arc). 



Figure 6.- Theoretical effect of change in thickness distribution on the acoustic pressure signature. 

Thickness noise only. (Note the change of scale.) 
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(b) Helicopter speed, 170 knots. 


Figure 7,- Comparison with experiments. 
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Figure Al.- The collapsing sphere g = 0 in the process of intersecting the 
moving body f(y,T) = 0 and forming the F-curve. The S -surface is the 
surface generated by all F-curves as the source time t varies from 
-°o to t for fixed observer position x and time t. 
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Wing velocity 



(a) Wing velocity. 



Observer 


(b) Wing position at time Tq. 


Figure A2. - Example of a condition causing singularity in the solution 

of the wave equation. 



Figure A3. - The S -surface corresponding to the observer time tq + Tq/c 
of the wing of figure A2 some of the F -curves are also shown. 


Wing velocity 



Figure A4,- Time history of the motion of the wing of figure A2 causing 
the S -surface to intersect itself. 



Figure A5.- S -surface of the wing of figure A2 corresponding to the observer time 
ro + Tq/c. Some of the F-curves are also shown. 
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Figure A6.- Structure of the 2 -surface near point B of figure A5. 





(a) Point A of figure A3. (b) Point B of figure A3. 

Figure A7.- The F- and m -curves in various situations causing singularities 
in (F,t) - formulation of the solution of the wave equation. 
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g = T - t + r/c = 0 



Figure A8.- Diagram for determination of the relation between the E -surface 
and (r,r) -formulations. AC = c dr; AB = dm. 



Figure A9.- Diagram for determination of the relation between the elements 
of the areas of the 2 -surface and the body surface S. (The unit vectors 
r and m are tangent to F-curve and m -curve, respectively.) 
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